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ABSTRACT 

We present a method for determining the B-field around neutron stars based 
on observed kHz and viscous QPO frequencies used in combination with the best- 
fit optical depth and temperature of a Comptonization model. In the framework 
of the transition layer QPO model, we analyze magnetoacoustic wave formation 
in the layer between a neutron star surface and the inner edge of a Keplerian 
disk. We derive formulas for the magnetoacoustic wave frequencies for different 
regimes of radial transition layer oscillations. We demonstrate that our model 
can use the QPO as a new kind of probe to determine the magnetic field strengths 
for 4U 1728-42, GX 340+0, and Sco X-1 in the zone where the QPOs occur. Ob- 
servations indicate that the dependence of the viscous frequency on the Keplerian 
frequency is closely related to the inferred dependence of the magnetoacoustic 
wave frequency on the Keplerian frequency for a dipole magnetic field. The 
magnetoacoustic wave dependence is based on a single parameter, the magnetic 
moment of the star as estimated from the field strength in the transition layer. 
The best-fit magnetic moment parameter is about (0.5 — 1) x 10^^ G cm^ for all 
studied sources. From observational data, the magnetic fields within distances 
less 20 km from neutron star for all three sources are strongly constrained to be 
dipole fields with the strengths 10''^^ G on the neutron star surface. 
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Subject headings: accretion disks — stars: neutron — stars: individual (4U 
1728-34, Scorpius X-1; GX 340+0) — stars: magnetic fields — X-rays: bursts 

1. Introduction 

Strong magnetic fields of neutron stars, typically of order ~ 10^^ G (Bhattacharya, & 
Srinivasan 1995) are considered crucial to the radio pulsar activity but equilibrium argu- 
ments concerning NS magnetic field pressure and ram pressure balance lead to magnetic 
field strength estimates for heavily accreting NS sources (for example X-ray bursters) of 
^ 10^ G. In these sources no persistent pulsations have been found despite careful searches 
(Wood et al. 1991; Vaughan et al. 1994). Different scenarios of magnetic field evolution 
lead to different estimates of magnetic field strength and decay. Thus, it is desirable to find 
a direct way of measuring or estimating the neutron star's magnetic field in these binary 
systems. Any model used to make such estimates should have minimum assumptions about 
observed data and model relationships to produce accurate results. Our model, described 
herein, has only one parameter that can be adjusted for agreement, the magnetic moment. 
All other parameters are strongly constrained by observations. It develops a new probe of 
the magnetic field strength and geometry, based on radial magnetoacoustic oscillations. 

Titarchuk, Lapidus & Muslimov (1998), hereafter TLM98, considered the possibility 
of dynamical adjustments of a Keplerian disk to the innermost sub-Keplerian boundary 
conditions to explain most observed QPOs in bright low mass X-ray binaries (LMXBs). 
They concluded that an isothermal transition sub-Keplerian layer between the NS surface 
and its last Keplerian orbit forms as a result of this adjustment. The TLM98 model is 
general treatment applicable to both NS and black hole (BH) systems. The primary problem 
in both NS and BH systems is understanding how the fiow changes from pure Keplerian 
to sub-Keplerian as a radius decreases to small values. The TLM98 authors suggested 
that the discontinuities and abrupt transitions in their solution result from derivatives of 
quantities such as angular velocities (weak shocks). Subsequently, Titarchuk & Osherovich 
(1999 hereafter, T099) identified the transition layer (TL) model's corresponding TL viscous 
frequency uy and diffusion or break frequency z/f, for 4U 1728-34, and predicted values for vt, 
related to diffusion in the boundary layer. TLM98 argued that a shock should occur where 
the Keplerian disk adjusts to sub-Keplerian flow. They interpreted the low frequencies as 
acoustic oscillations and expected them to provide reliable estimates of radial velocities that 
would be very close to the acoustic velocities. 

Sco X-l's derived acoustic velocities (Titarchuk et al. 2001 hereafter TBGF) translated 
to a TL plasma temperature, assuming only thermal motion, of about 19 keV. However 
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photon spectral fitting by TBGF yielded a TL temperature of 2.8 keV. TBGF attributed 
this discrepancy to a dominant contribution to the plasma matter energy density from the 
magnetic field. The viscous QPOs provide information identifiable with a magnetoacoustic 
(MA) oscillation effect that can be modeled. 

In this Letter, we develop a general formahsm for MA radial oscillations in the TL. In 
§2, we describe the data sources used to test our MA oscillation model. The formulation 
of the problem and derivation of MA frequencies are described in §3. We show the results 
of the MA oscillation interpretation for the T099 QPO viscous frequencies and the best- 
fit normalization parameters of magnetic field strengths for a number of QPO sources are 
present in §4. Our summary and conclusions follow in §5. 

2. Observations 

The model can only be applied to QPOs in which the full TLM has been worked out 
and in which the QPOs identified with oscillations of a bounded medium have been seen. All 
observations^ referenced herein were made with the Rossi X-ray Timing Explorer's (RXTE). 
Sco X-1 data were obtained during August 3 and 22, 1997, February 27 and 28, 1998, and 
June 10-13, 1999. Data were extracted from research archives for GX 340-1-0 and 4U 1728-34 
for observations on November 14, 1998 and February 28, 1999, respectively. In addition, we 
use the results from observations of 4U 1728-32 between February 15-March 1, 1996 by Ford 
& van der Khs (1998) 

3. Magnetoacoustic Oscillations in the Transition Layer 

In the references cited above, the treatment of a magnetic field in a fluid disk has 
specifled neither the multipole order that characterizes the held nor the strength of the held. 

We derive the frequency of the QPO associated with MA oscillations and the correlation 
of the MA frequency with the Keplerian frequency u^. The MA frequency is derived as the 
eigenfrequency of the boundary- value problem resulting from a MHD treatment of the inter- 
action of the disk with the magnetic fleld. The problem is solved for two limiting boundary 
conditions encompassing reahstic possibilities. The solution yields a velocity identified as a 
mixture of the sound speed and the Alfven velocity, becoming either at the appropriate lim- 



"'^This research has made use of data obtained through the High Energy Astrophysics Science Archive 
Research Center Online Service, provided by the NASA/Goddard Space Flight Center. 
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its. Our treatment does not specify how the eigenfrequency is excited or damped. However, 
it makes clear that the QPO is a readily stimulated resonant frequency. 

In order to derive the MA oscillation frequencies in the isothermal transition layer, we 
consider layer perturbations such that the density p — po + pi, the velocity v = Vq + Vi and 
B = Bo + Bi. By combining the continuity equation, equation of motion, and an ideal gas 
equation, we obtain an equation for MA oscillations (Alfven 1942) 

- s2V(V • vi) + VA X V X [V X (vi X va)] = F(vo, dvo/dr, Bq, po) (1) 

where = {dP/ dp) is the square of the sound velocity (assumed constant with respect to 
radius) and Va = Bo(47rpo)"''^^^ is the Alfven velocity. To determine the eigenfrequencies of 
MA oscillations in the TL, we set equation (1) to zero and introduce a cylindrical coordinate 
system in the TL with a radial axis in the disk plane, z-axis perpendicular to the plane, 
and azimuthal axis along the circle in the disk plane [see Fig. 1 in Titarchuk, Osherovich & 
Kuznetsov (1999), hereafter TOK, for the transition layer geometry]. 

The magnetic field acts along a direction in space that is the same throughout the TL, 
intersecting the disk everywhere at the same angle 5. Practical ranges for 5 are from a few 
to ~ 15 degrees (Titarchuk & Osherovich 2001). The model used for all three sources (see 
§4 ) assumes that 6 = exactly, which reinstates axial symmetry and makes the problem 
mathematically tractable. The simplified problem is spatially one dimensional in the radial 
direction, and oscillations represented by the eigenfunctions are radial fluid movements. 

We consider the radial oscillations in the TL for which we should assume that Bq = B^e^, 
i.e., the mean magnetic field strength Bq is perpendicular to the disk and Vi = Vie,. is along 
the radial direction in the disk. We wish to explore how solutions for fiow in the TL are 
affected by boundary conditions and B field geometry, for acoustic and magnetic oscillations. 
Under these assumptions, equation (1) can be rewritten as: 

d'^vi „2 ^ A ^^^1 \ y ^ ("^ dry aVi \ _q 
dt"^ dr \r dr J dr \r dr J 

Combining this equation with the two boundary conditions of the TL, the neutron star 
surface r = it!j„ and the outer boundary, r = Routi defined by the Keplerian QPO frequency, 
allows us to analyze the eigenfrequencies of the MA oscillations in the TL under both the stiff 
vi{Rin) — vi{Rout) — and the free dvi/dr{Rin) = dvi/dr{Rout) — boundary conditions. 

In the two extreme cases, either va <^ s or v\ ^ s equation (2) can be examined 
analytically under the appropriate boundary conditions. Introducing a new function y = rvi 
in the case oi va <^ s and, y = tvaVi, in the case of va ^ s, we can write equation (2) in 
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the form 

where x is a new variable, a — and A — As = s"^ for the pure acoustic case, and a > 6 
and A — A^ — B'^r^/Anp for the magnetic case. The index a is related to the multipole 
magnetic field, v\ — A^/r"" {a — 6,8,10 are for the dipole, quadrupole, and octopole, 
respectively). is the magnetic field strength at a radius < r* < Rout, ^-nd p is the TL 
mean density. 

First, we need to find a solution of the eigenvalue problem for equation (3) under the 
stiff boundary conditions by obtaining a particular nontrivial solution using the separation 
variables method, such that 

y{x,t) ^ X^{x)T^{t) (4) 
where Xx{x) is an eigenf unction satisfying 

(a-2)/2 

X" + \ , X^X = (5) 

4:A ^ ' 

and boundary conditions X{xi^ — X{xout) — 0. The solution of these equations can be 
expressed through the Bessel Z functions 



2 ^ .2.(a+2)/4 



(6) 



{a + 2) AV2 

X\ is a linear superposition of the Bessel ^2/(0+2) and J_2/(a+2) functions: 

Xx^x'/^[J2/(a+2){zx) + CJ_2/{a+2){zx)] (7) 

where 

' (a + 2)AV2^ ■ i«i 

Wood et al. (2001) demonstrated that an asymptotic form of the Bessel function for zx^l 
can be used to find all eigenvalues including the first one: 

/ 2 \ 

^2/(a+2) (zx) = cos[zx - (« + 6)n/4:{a + 2)] (9) 

\7rzx/ 

and 

/ 2 

J-2/ia+2)(^x) = cos[zx -{a- 2)77 /4{a + 2)]. (10) 

\^7lZx J 

Using these asymptotic forms, the boundary conditions, and equations (7 & 8), the 
eigenvalues are determined by 

Q,0^{Zout - </?l) ^0%{Zin - <P2) " COs{Zout " ^2) COs{Zin - <P\) = 0, (11) 
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where cpi — {a — 2)7r/4(Q; + 2) and (p2 — {a + 6)71/4(0; + 2). Equation (11) can be reduced 
to sm{zout ~ Zin) — or (3 — Zout ~ ^in — ^TJ"; ^ = 1) 2, 3..., with a solution 



2 „(«+2)/2_ (a+2)/2- (12) 
out ' in 



The time dependence, T(t) is found from the equation T" + A^T = 0. Thus, T{t) = sin(A„i + 
0o) = sin(27rz/„t + 0o) where 0o is an initial phase. For the acoustic case (A = Ag = s and 
a = 0), the main eigenvalue Ai = TT/{rout — ^in) = t^I L and I's = 1^1 = Ai/27r = s/2L, where 
L is the TL radial size. 

For the magnetic case (a > 6) 

(« + 2)5, / r, \ '^Z' (a + 2)vj^{rout) „n 



4(471 py/^[r out - rin{rin/routY^'^] \routJ ^[rout - rin{rin/routY''^ \ 

For the free boundaries, the boundary conditions for y are 

dy , {a-2) y 

— H = aX. x = Xin and a; = Xout, (14) 

from which, the transcendental equation for determining the eigenvalue is 

^^"^ - - (« + 2){r?/5V(r/ - 1)2 + [(a - 2)/{a + 2)]^} ' ^ ^ 

where p = 2;o„t - Zin and 77 = Zout/zin = (rout/rj„) ("+2^/^ 

Analysis of equation (15) shows that for the values most relevant to observations (77^2), 

/3M«^^ + -[(7r/2)V(^-lf + l/4] (16) 

for a = 6 and (3^ ^ {1.5/[1 + 1.5r]/{r] - 1)2]}V2 for a = 0. 

The free-boundary conditions will therefore introduce additional factors appearing in 
formulas for Ui between 0.5 and 1 for the magnetic case and I/tt for the acoustic case. From 
this, we conclude that the free-boundary conditions can decrease substantially the magnetic 
oscillation frequency and must decrease the acoustic frequency by at least a factor of 3 (see 
Eqs. 17-18 below). 
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4. Magnetic Field Strength Determination and the Magnetoacoustic 

Oscillation Frequency 

Wc can construct an approximate formula for the MA frequency i'ma using the asymp- 
totic forms for the two extreme (acoustic and magnetic) cases 

l^MA ~ {l^'s + ^mY'^ = {s^/^{r^t - Tin? + [(« + I vI{t ^t) / [t - Ur^iViJ T^tY/^ff/^ 

(17) 

The approximate relation in equation (17) is exact if the Alfven velocity is assumed to be 
constant through the TL. The difference, by a factor (a + 2)^/4, in coefficients for and v\ 
is due to different behaviors of s and vx as a function of r (s = const and va oc r~"/^). 

Under the free boundary, equation (17) is modified to 

yuA ^ m/^fs^/Krout - Un? + (/3M/7r)^[(a + 2)/^fvl{Tout)/[T^t - rin{rin/routr"YY'^ 

(18) 

and provides the dependence of vma on i/k- This dependence is due to the explicit dependence 
of Tout on the Keplerian frequency uk: 

r^t = (GM/27r4)V3, (19) 

where G is the gravitational constant and M is the NS mass. To determine s, wc define 
s = [kT /rripY/'^ as the sound speed for the isothermal TL where kT is the best-fit value 
of the Comptonization temperature. To determine vm = B{rout)/ {'^'^PoY^'^ ■, we express 
p — rnip/ {(TtL) with r as the best-fit optical depth parameter. These are not adjustable 
parameters for determining magnetic field strengths. 

Thus, we obtain Vma a function of only one adjustable parameter, the magnetic- 
field strength at a TL radius r*. We can replace this parameter by the multipole moment 
B{r^,)r"^^. In the dipole case a — 6, the multipole moment is the dipole moment (in G 
cm^). The procedure for determining the magnetic field strength in the TL may reahzed if 
there is an observable dependence between the viscous and Keplerian frequencies. 

The spectra of 4U 1728-34, GX 340+0, and Sco X-1 were fitted by the Comptonization 
model of Sunyaev & Titarchuk (1980) from 8-25 keV, to obtain a plasma temperature (kT) 
and Thompson optical depth (r). Figure 1 shows the results from the model and the fitted 
data for 4U 1728-34, GX 340+0, and Sco X-1 (Ford & van der Klis 1998). 

We found that a dipole geometry (a = 6) gives a good fit and higher poles are excluded 
by observed data. It is also evident from equations (17-19) that the MA frequencies are well 
described by a power law dependence on when vm > t^s, i-e., i'ma ~ t^m oc u^'^^^^^. The 
power law index in the dipole case is (a + 2)/3 = 8/3. The MA frequency dependence on 
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vk is too steep for the quadrupole and the octopole cases [where {a + 2)/3 = 10/3 and 4, 
respectively] to be supported by the data. The fits of our MA model to the data provide 
strong constraints for the boundary condition types where the stiff conditions are ruled out 
for 1728-34 and the free conditions for GX 340+0. The fit qualities are practically insensitive 
for M = 1.2 ± O.2M0. The best-fit mass parameter is around 1.2. 

The best-fit parameters for 4U 1728-34 (Compton values kT = 8.6±2 keV, r = 5.8±0.5; 
x'^/dof = 19.8/35) are = (0.86 ±0.03) x 10^ G, and r, = 1.69 x 10^ cm, corresponding to 
a magnetic moment /i = 4.2 x 10^"^ G cm^. We fixed the NS radius, rj„ to be (as a canonical 
value) three Schwarszchild radii 3 r^. An extrapolation of the magnetic field from towards 
the NS radius gives us B^s = 4.5 x 10® G and B^s = 1.3 x 10^ G for the dipole and 
octopole fields, respectively. 

The best-fit parameters for GX 340+0 (Compton values kT = 2.87 ± 0.02 keV, r = 

11.2 ± 0.2; x^/dof = 49.8/35) are B^ = (0.5 ± 0.05) x 10® G and = 2.7 x 10® cm, 
corresponding to a magnetic moment = 1 x 10^^ G cm^. An extrapolation of the magnetic 
field from r* towards the NS radius rj„ yields B^g — 8-3 x 10^ G and B^s — 5.3 x 10^ G for 
the dipole and octopole fields, respectively. 

The best-fit parameters for Sco X-1 (Compton values kT = 2.87 ± 0.03 keV, r = 

11.3 ± 0.3; x^/dof = 9.9/33) are 5, = (1. ± 0.05) x 10^ G, and n = 2.12 x 10^ cm, 
corresponding to a magnetic moment /j, — 9.5 x 10^^ G cm^ in the TL. An extrapolation 
of the magnetic field from r* towards the NS radius ri„ gives us B^s — 8-1 x 10^ G and 
Bns — 3.3 X 10^ G for the dipole and octopole fields, respectively. 

These derived magnetic-field strengths are for the TL only and can be smaller than the 
true magnetic field outside of the layer because of its finite conductivity. The NS magnetic- 
field strengths can be increased by factors of 1.7 and 2.5 for the dipole and octopole extrap- 
olation respectively if we assume that the NS radius is around 2.5 rs instead of 3 rs- 

The magnetic field strength can be roughly checked using an equipartition estimate for 
which S^/Stt = nkT. For example, Sco X-l's TL density p = 2.7 x 10~^ g cm~^ is obtained 
using the best-fit optical depth r = 11.3. The TL size, L — 10^ cm, is obtained using the 
observed Keplerian frequency and best-fit temperature kT — 2.9 keV. Thus, we obtain a 
magnetic field, S = 1.3 x 10® G, at radius r = 2.3 x 10® cm. 
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5. Summary and Conclusions 

Since the discovery of neuron stars in the 1960s magnetic fields have been estimated 
almost exclusively from rotational effects, i.e. from periodic pulses with their periods p and 
derivatives p [e.g. Bhattacharya & Srinivasan (1995)]. The new idea here is that, when the 
field configuration is nearly axis-symmetric, there may arise circumstances in which the field 
is estimated from radial oscillations (QPOs) of plasma interacting with magnetosphere. 

The model preconditions are: (1) the MA QPO must be detected observationally, (2) a 
Keplerian QPO frequency must be observed and (3) the temperature and the optical depth 
of the QPO emission region must be obtained from the photon spectra. Prom our model 
we derived tightly constrained dipole magnetic fields of the TL for the bright LMXBs 4U 
1728-34, GX 340+0, and Sco X-1. The observed low frequency correlation of the kHz QPOs 
strongly confine the theoretical dependence of uma on z/k, which is mostly determined by 
the magnetic frequency dependence on z/k (see dashed lines in Figure 1) and very sensitive 
to the pole order (dipole, quadrupole, etc.). The correction to acoustic oscillations [the first 
term in formulas (17-18)] only slightly changes the theoretical shape of the low frequencies 

of Z/R. 

We re-fit the stellar magnetic dipole moment in each source to derive a range of values 
from (0.5 — 1) x 10^'^ G cm"^. The field estimate arrived at this manner is close to that may 
be evolutionary end points of LMXBs, namely milisecond pulsars, which have derived values 
(from v^) of ~ 10^ G (Bhattacharya 1995). 

We wish to thank Sergey Kuznetsov, Eric Ford, Peter Jonker and Michael van der Klis 
for the Sco X-1, 4U 1728-34 and GX 340+0 data files. We are truly grateful to the referee, 
J. Imamura for his thorough analysis of the presented manuscript. 
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Fig. 1. — MA frequency vs Keplerian frequency for 4U 1728-34 (black squares represent 
QPO frequencies from Ford & van der Klis 1998, blue, from TOK and magenta, from Jonker 
et al. 2000). 1-a errors bars indicate statistical fitting errors (A^^ = 1). The value for Q (= 
Z//FWHM) approaches 1.0 for viscous QPOs corresponding to Keplerian frequencies of 320- 
360 Hz (Ford & van der Klis 1998), giving formal frequency errors significantly larger than 
the plotted error bars. The solid lines are the theoretical best-fits for the MA frequencies. 
The dashed lines are for the magnetic frequencies only (see text). 



